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Abstract 

In this work we present some classes of models whose the corre- 
sponding two coupled first-order nonlinear equations can be put into 
a linear form, and consequently be solved completely. In these cases 
the so called trial orbit method is completely unnecessary. We recall 
that some physically important models as, for instance, the problem 
of tiling a plane with a network of defects and polymer properties are 
in this class of models. 
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A rapid look at the history of physics is enough to lead anyone to con- 
clude that, fortunately, the most part of the natural physical systems can 
be studied by using linear differential equations, with their good properties 
like the superposition principle. Notwithstanding, there are some classes of 
important systems with are intrinsically nonlinear and, nowadays, there is a 
growing interest in dealing with such systems pQ - [TH]. Unfortunately, as a 
consequence of the nonlinearity, in general we lose the capability of getting 
the complete solutions. In this work we show that for those systems in 1+1 
dimensions, whose the second-order differential equations can be reduced 
to the solution of corresponding first-order equations, the so called Bolgo- 
mol 'nyi-Prasad-Sommerfield (BPS) topological solitons [T7], one can obtain 
a differential equation relating the two coupled fields which, once solved, 
leads to the general orbit connecting the vacua of the model. In fact, the 
"trial and error " methods historically arose as a consequence of the intrin- 
sic difficulty of getting general methods of solution for nonlinear differential 
equations. About two decades ago, Rajaraman [18] introduced an approach 
of this nature for the treatment of coupled relativistic scalar field theories in 
1+1 dimensions. His procedure was model independent and could be used for 
the search of solutions in arbitrary coupled scalar models in 1+1 dimensions. 
However, the method is limited in terms of the generality of the solutions 
obtained and is convenient and profitable only for some particular, but im- 
portant, cases [15]. Some years later, Bazeia and collaborators [20J applied 
the approach developed by Rajaraman to special cases where the solution of 
the nonlinear second-order differential equations are equivalent to the solu- 
tion of corresponding first-order nonlinear coupled differential equations. By 
the way, Bazeia and collaborators wisely applied their solution to a variety 
of natural systems, since polymers up to domain walls. In this work we are 
going to present a procedure which is absolutely general when applied to 
systems like those described in [20], namely the BPS topological solutions. 
Furthermore, we are going also to show that many of the systems studied 
in [20]- [25] can be mapped into a first-order linear differential equation and, 
as a consequence, can be solved in order to get the general solution of the 
system. After that, we trace some comments about the consequences coming 
from these general solutions. 

In order to deal with the problem, following the usual procedure to get 
BPS [T7] solutions for nonlinear systems, one can particularize the form of 
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the Lagrangian density 



(1) 



by imposing that the potential must be written in terms of a superpotential 
like 



1 ( dW(<f>, X ) 



2 + ifdw^, x )\ 2 
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The energy of the so-called BPS states can be calculated straightfor- 
wardly, giving 
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which lead us to 
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and finally to 
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where <pi and Xi are the are the i — th vacuum state of the model [22] . 

In this case, one can easily see that solutions with minimal energy of 
the second-order differential equation for the static solutions in 1+1 dimen- 
sions, can be solved through the corresponding first-order coupled nonlinear 
equations 

% = w*(M,^ = w x (<i>,x), (6) 

where = ^ and W x = Here, it is important to remark that the 
BPS solutions settle into vacuum states asymptotically In other words, the 
vaccum states act as implicit boundary conditions of the BPS equations. 

Now, instead of applying the usual trial-orbit approach [2H]-[25], we note 
that it is possible to write the following equation 



dx 



dx 



(7) 
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where the spatial differential element is a kind of invariant. So, one obtains 
that 

d4 = W ± . , 

d X W x 1 ; 

This last equation is, in general, a nonlinear differential equation relating 
the scalar fields of the model. Now, if one is able to solve it completely, 
the function <p (x) can be used to eliminate one of the fields, so rendering 
the equations ([6]) uncoupled. Finally, this uncoupled first-order nonlinear 
equation can be solved in general, even if numerically. 

From now on, we choose a particular model which can be used for mod- 
elling a number of systems [22], in order to exemplify the method in a concrete 
situation. In fact we will show that for this situation, the equation (jSJ) can 
be mapped into a linear differential equation, from which it is possible to 
obtain the general solution. In this case the superpotential is written as 

iy(</>, X ) = -A0+^ 3 + /x0x 2 , (9) 
and the equation (jSJ) looks like 

d4> A(0 2 -1) + ^ X 2 



dx 2 n (j) x 



(10) 



At this point one can verify that, performing the transformation 2 = 
p + 1. The above equation can be written as 

dp X 

-i P = X> 11 

dx 

a typical inhomogeneous linear differential equation. It is interesting to ob- 
serve that its particular solution, corresponds to the result usually presented 
in the literature [22]. The general solution is easily obtained, giving 

p ( X ) = 2 - 1 = c x$ ~ tAt X 2 ' ( 12 ) 

A — Z p 

for A 7^ 2 p, and 

2 -l = X 2 [ln(x) + Cl ], (13) 

for the A = 2 p case, and Cq and C\ are arbitrary integration constants. It 
is interesting to note the this last particular situation was not taken into 
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consideration in the literature up to our knowledge. From now on, we sub- 
stitute these solutions in one of the equations (j5D, and solve it, so obtaining 
a generalized solution for the system. In general it is not possible to solve \ 
in terms of <j) from the above solutions, but the contrary is always granted. 
Here we will substitute (j) (%) in the equation for the field x, obtaining 



ax \ 



1 + CoX ^~\X^2jl) x2 ' ^ 2 ^' 



and 



^ = ±2fj, X yJl + X 2 [ln(x) + ci], (A = 2/i). (15) 

In general we can not have an explicit solution for the above equations. 
However one can verify numerically that the solutions are always of the same 
classes. Notwithstanding, some classes of solutions can be written in closed 
explicit forms. First of all, we should treat the system when cq = 0, because 
in this situation we can solve analytically the system for any value of A, apart 
from the case A = 2 \x. In this situation we get 

2 g 2/i (z-zo) 2 e 4 K x ~ x o) 

^+ = l _ ce i^(x-x ) ' X- ( X ' = c _ e 4fi(x-x ) ' 

with c = — A _ M 2 - . For this choice of the parameters, the solution always 
vanishes at the boundary (x — > ±oo). As a consequence, the corresponding 
kink solution for the field 0, will be given by 

ce 4n(x-x ) _|_ j c _|_ e 4fi(x-x ) 

^ ^ = ± ce 4„(s-*o)-l ' ^ ^ = ± c _ e 4^-x ) » 

which are essentially equivalent to those solutions appearing in [22], given in 
terms of tanh(x). Let us now discuss below two particular cases (c ^ 0) 
where the integration can be performed analytically up to the end. Let us 
first consider the case were A = fi, which has as solutions 

^ e 2fj,(x-x ) 4 g2/i(a;-xo) 

X+ ^ ~~ [ Cq e M*sa) - I] 2 _4 e 4M(s-a:o)' X ~ ^ ~~ [ e 2^(x-x ) _ Cq ] 2 _ 4' 

(18) 

where we must impose that Co < — 2 in both solutions, in order to avoid 
singularities of the field as can be easily verified. Furthermore, both solutions 
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vanishes when x — > ±00, provided that Cq 7^ — 2. On the other hand the 
corresponding solutions for the field <fi (x) are given by 

(eg - 4) e 4 ^ x - x °) - 1 , . . 4 - cl + e 4 ^- x °) 

W = 7 — m ; — wi — r~n 7; <P- ( x ) ~ 



[Co e 2 ^ x ~ x o) — l] 2 — 4 e 4v(x-xo) ' [g2|*(a:-xo) _ C() j 2 _ 4 

(19) 

Here the first bonus coming from the complete exact solution of the equa- 
tion ([6]) comes when we deal with the special case with cq = —2. It is re- 
markable that for this precise value of the arbitrary integration constant, 
an absolutely unexpected kink solution do appears. In fact, it could never 
be obtained from the usually used solution, where cq = necessarily In 
this special solution, the field x is a kink with the following asymptotic lim- 
its: x+ ( — 00 ) = and x+ (°°) = lj an d 4>+ ( — 00 ) = ~~ 1 an d 4>+ (00) = 0, 
and correspondingly X- (~ 00 ) = 1 an d X- (°°) = 0, and (—00) = and 
0_ (00) = 1, as it can be seen from an example of a typical profile of this kink 
in the Fig. 1 . Below we present a plot of this kink, which we are going to call 
type B kink, in contrast with the other cases where the field x does not have 
a kink profile, which we call type A kink (see Fig. 2) . An interesting obser- 
vation is that the choice cq = —2, is precisely the one which makes the right 
hand side of equation (TI4"]) simply proportional to x 1 1 — xl = Cx (1 ~ x)> 
where ( is the sign function defined as ( = (1 — x) 1 11 — x\ ■ It takes values 
±1 with £ = +1 being selected by boundary conditions < x ^ 1 f° r the 
solutions appearing in (18) and, in this situation, the equation is much easier 
to solve. In fact, by performing the translation x — P + \, we recover a BPS 
superpotential for the "A0 4 " model, — (/3 2 — 1/4). A similar situation will 
happens with the next example. 

As the third particular case, we consider the situation where A = 4//. 
Now, the exact solutions look like 



x+ 0) 



x- 0) 



2 e 2 K x ~ x o) 

"' 1 

■ Cq e 8 ^ -1 ' ) 

(20) 



i 2 

I e 4[t(x-x ) _|_ ]_ — 4 Cq e 8 ^( x_x °) 



^_ e 2ti(x-x ) 



1 + 2 e 4 ^"")] 2 - 16 c ' 



which have the same asymptotic behavior as that presented in the previous 
cases for the type A kinks. In other words, provided that Cq 7^ 1/16, only 
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the field <p will be a kink. Afterwards, as in the previous case, if one wish to 
avoid intermediary singularities, one must impose that Cq < Now, the 4> 
solutions will be written as 

4+ (16 c - 1) e 8 ^^ ) 
[2 + e 4 ^ x - x °)] 2 - 16 c e 8 ^" 10 ) ' 

(21) 

16 c + 4e 8 ^- x °) - 1 
[1 + 2 e 4 ^^o)] 2 - 16 c ' 

Once more, the particular choice of the integration parameter Co = jg, 
generates a type B kink, with the asymptotic behavior given by: x+ (~ 00 ) — 
and x+ (°°) — — 2, and (fi + (— oo) = 1 and (f> + (oo) = 0, and correspondingly 
X- (— oo) = —2 and (oo) = 0, and 0_ (— oo) = and 0_ (oo) = 1. 

It is interesting to calculate the energy of these two species of solitonic 
configurations. For this we use the superpotential ([9]) and substitute it in 
the equation ([5]), and observe that the type A kinks have an energy given 
by E^ PS — I A and in the two cases considered above (A = /i and A = 4 jj) 
we obtain E§ ps = | A. One could interpret these solutions as representing 
two kinds of torsion in a chain, represented through an orthogonal set of 
coordinates 4> and %. So that, in the plane {(fi,x), the type A kink corresponds 
to a complete torsion going from (—1, 0) to (0, 0) and the type B corresponds 
to a half torsion, where the system goes from (—1,0) to (0, 1), in the case 
where (A = fi) for instance. 

In what follows, we will study a more general model, contemplating a 
number of particular cases which have been studied in the literature, includ- 
ing the previous and some other new ones. For this, we begin by defining 
the superpotential 

W^, X ) = ^ N X 2 + G( ( p), (22) 
which lead us to the following set of equations: 

f x = d - £ ^ + l N ^ 1> ^ 1 = "*"* < 23) 

So, the corresponding equation for the dependence of the field as a function 
of the field Xi is given by 
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Now, performing the transformation a = <p 2 we get 



da = N 1 



(25) 



where G t / ) (a) = | cr=0 2 - Obviously, there are no arbitrary solutions for 

the above equation, but for that ones with exact solution we can get the 
corresponding exact two-field solitons. For instance, let us treat the special 
case where 



dG{<j>) 



2 (an + ai o + a 2 a 2 ) n=i 
— -a 2 . 



(26) 



The solution will be given by a combination of Bessel functions which, once 
substituted in the equation for the field x, lead us to a hardly exactly solvable 
equation, beyond some singularities which appear in the solution. So, we still 
here continue to work with the simpler linear case of this equation, where 
a2 = 0, which furthermore permits us to write arbitrary solutions given by 



v KX) = n( x 

d\ 2 (/i — CLi) 



2 ax 

+ ci x M , 



with G (0) given by 



G{<P) 



a 



+ 



N (N + 2) 



(27) 



(28) 



leaving us with the following potential 



V(o.\) = y;r x ^a\,y ■ 2r/,.V / /,r\ 2 • af t ■ // 2 \ 2 (l. 



+2a (ai0 2 + Nfi X 2 ) , 



+ N\ 2 ) + 

(29) 



with ci being the integration arbitrary constant, and a and ai are constants 
which characterize the physical system. From above, it is easy to conclude 
that 



= ± 



\ 



Oo 



Nnx 2 



2 ax 

+ ci x~, 



2(ji - a x ) 

and, consequently we are left to solve the following equation 

Nfi X 2 



— = ±fi 

dx 



a 



ai 2 (fj, — ai) 



2 ax 



(30) 



(31) 
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At this point it is important to remark that many models appearing in 
the literature can be cast as particular cases from the above general one. For 
instance if we take N — 1, we recover the models I, II and III of [22], and 
model I of [20]. The case where N = 2 is equivalent to the model II in [20] 
and the model considered in [23]. 

As a final comment we should say that one can even make a bit general- 
ization of the above exactly solved two fields models. This could be done by 
starting with the superpotential 



W NM (0, X) 



(32) 



with G (<f)) being the same appearing previously in the text. After manipu- 
lations similar to that one done above, we end with the equation 



da{x) 
dx 



'2a£ 



o (X) X (1 ~ M) + ( 



2N 



X ■ 



(33) 



where a = 2 + (^j . Solving the above equation for arbitrary M, one obtains 
that 



o"m (x) 



2aiX< 2 - A/ ) 
^(Af-2) 



2aiX 



{2- MY 



yU (2 — M) 



ci + 



1 



M(M-2) 



aix 



fi (2 - M) 



2m=2 n x 2 r 



(2-M) \ \-M=2, 



2 

M-2' 



(34) 



where c\ is the arbitrary integration constant, and T (a,z) = J z °° t^^e^dt, 
is the incomplete gamma function. 

Obviously, the case studied earlier in this work is obtained from the above 
when one chooses M — 2. On the other hand, we can get simpler solutions 
for other particular values of the parameter M as, for instance M — 4, whose 
solution can be written as 



^4 (X) = jX + e 



, N ai 
ci + — — Ei 

4yU 



(35) 



where i?i (2;) 



X^z i s ^ e exponential integral function. It can 



be seen from Fig. 3 that, apart from a small region close to the origin, it is 
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asymptotically similar to that of the case with M — 2, which was discussed 
in some detail above in the text. This expression does not have any kind of 
singularity and approaches to zero when the field x does the same. Notwith- 
standing, the last part of the analysis of the kinks needs to be done through 
evaluation of the equation 

AT 

(36) 

which is not easy to be done analytically, so that one needs to make use of 
numerical techniques. We intend to perform this analysis in a future work, 
looking for new interesting features. 



ax 



(M-l) 



(±) 



N 



— ) + <?M (x) 
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Figure 1: Typical type B kink profile (for X — fj, ). The dotted line corre- 
sponds to the field x+( x ) an d the solid line to the field <f> + (x). Both were 
calculated for c = — 2. 




Figure 2: Typical type A kink profile (for A = /i ). The dotted line corre- 
sponds to the field x+( x ) an d the ^o\id line to the field <f>+(x). Both were 
calculated for c = — 3. 
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